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We report an improved measurement of D + — > K%K + and Df — » Kgir + branching ratios using 
605 fb^ 1 of data collected with the Belle detector at the KEKB asymmetric-energy e + e~ col- 
lider. The measured branching ratios with respect to the Cabibbo-favored modes are B(D + — > 
K%K + )/B{D + -> K° s tt+) = 0.1899±0.0011±0.0022 and B{Dt -> K° s ir + )/B(Dt -> K° S K+) = 
0.0803±0.0024±0.0019 where the first uncertainties are statistical and the second are systematic. 
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Decays of charmed mesons play an important role 
in understanding the sources of SU(3) flavor symme- 
try breaking Such a breaking can originate from 
strong final-state interactions or interference between 
amplitudes with the same final state. In particular, 
D+ -> K°K+ and D+ -> K°tt+ are Cabibbo- 
suppressed (CS) decays that involve color-favored tree, 
annihilation, and penguin diagrams. For D + decays, the 

branching ratio B(D+ -> K (> K + )/B{D+ -> K°ir + ) de- 
viates from the naive tan 2 9c expectation [3j , due to the 
destructive interference between color-favored and color- 
suppressed amplitudes in D + — > K tt + [4|. However, 
converting experimental measurements of D decays that 
include K s branching ratios to those involving K° or K 
is not straightforward due to the interference between the 
doubly Cabibbo-suppressed (DCS) and Cabibbo-favored 
(CF) decay modes where the interference phase is un- 
known d, @|. In D+ decays to K° K+ and K°ir + final 
states, the ratio of the CS decay to the corresponding 
CF decay may be larger than tan 2 6c, since the tree dia- 
gram for Df — > K K + is CF but color-suppressed. Pre- 
cise measurements of branching ratios for CS and CF 
charmed meson decay modes can thus improve the un- 
derstanding of the underlying dynamics of these decays. 
In this paper, we report improved measurements of the 
D + — > K S K + and — > K s tt + branching ratios with 
respect to the corresponding CF modes, D + — > Kgir + 
and Df — » K S K + , respectively. 

The results are based on a data sample of 605 fb _1 
recorded at the T(4S) resonance with the Belle detec- 
tor at the KEKB asymmetric-energy e + e~ collider Q. 
An additional data sample with about 10% of this inte- 
grated luminosity recorded 60 MeV below the T(4S) was 
used for the optimization of the selection criteria (off- 
resonance sample). The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 



(ACC), a barrel- like arrangement of timc-of- flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic 
field. An iron flux return located outside the coil is in- 
strumented to detect K\ mesons and to identify muons. 
The detector is described in detail elsewhere Q. 

We require that the charged tracks originate from the 
vicinity of the interaction point (IP) with the impact 
parameters in the beam direction (z-axis) and perpen- 
dicular to it of less than 4 cm and 2 cm, respectively. 
All charged tracks except those originating from Kg de- 
cays are required to have at least two associated hits 
in the SVD, both in the z and radial directions, to as- 
sure good spatial resolution on the D mesons' decay ver- 
tices. Charged tracks are identified as pions or kaons by 
requiring the ratio of particle identification likelihoods, 
£k/{Ck + At)) constructed using information from the 
CDC, TOF, and ACC, be larger or smaller than 0.6, re- 
spectively. For both kaons and pions, the efficiencies and 
misidentification probabilities are 86% and 10%, respec- 
tively. 

Pairs of oppositely charged tracks that have an invari- 
ant mass within 30 MeV/c 2 of the nominal Kg mass are 
used to reconstruct K s — > 7r + 7r~ decays. The distance 
of the closest approach of the candidate charged tracks 
to the IP in the plane perpendicular to the z axis is re- 
quired to be larger than 0.02 cm for high-momentum (> 
1.5 GeV/c) Kg candidates and 0.03 cm for those with 
momentum less than 1.5 GeV/c. The 7r + 7r~ vertex is re- 
quired to be displaced from the IP by a minimum trans- 
verse distance of 0.22 cm for high-momentum candidates 
and 0.08 cm for the remaining candidates. The mismatch 
in the z direction at the K s vertex point for the tt + tt~ 
tracks must be less than 2.4 cm for high-momentum can- 
didates and 1.8 cm for the remaining candidates. The di- 
rection of the pion pair momentum must also agree with 
the direction of the vertex point from the IP to within 
0.03 rad for high-momentum candidates and to within 
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0.1 rad for the remaining candidates. 

These two sets of criteria in two different momentum 
ranges are implemented to maximize Ms/V-^s +Mb, 
where Ms and Mb are the number of signal K s 's and 
the number of combinatorial background events, respec- 
tively. Finally, the ir + ir~ pair forming a K s candidate is 
required to have an invariant mass within ±9 MeV/c 2 of 
the nominal Kg mass [3J. 

D + and D+ candidates are reconstructed using a Kg 
candidate and either a charged pion or kaon candidate. 
The decay vertex is formed by fitting the Kg and the 
track to a common vertex and requiring a confidence level 
greater than 0.1%. In order to remove peaking back- 
grounds from the D7. — > 7r + 7r + 7r~ and K + tt + tt~ decay 

modes, we compute the reduced x 2 of the vertex assum- 
ing that two pions from the K s and the charm daughter 
track arise from a single vertex. We require the reduced 
X 2 to be greater than 10. 

To remove D + and Dj mesons produced in B meson 
decays, we require the charmed meson momentum cal- 
culated in the center-of-mass frame to be greater than 
2.6 GeV/c. At this stage, reconstruction efficiencies are 
16.6% for the D+ and 18.0% for the D+ in the K%K+ 
final state, and 20.6% for the D+ and 22.4% for the D+ 
in the K { gir + final state. 

Highly asymmetrical Kgh + pairs that have invari- 
ant mass close to the Df s ^ mass region are more likely 
to be background than signal. The asymmetry, A = 
\Pk° s ~ Ph+ \/\Pk° + Ph+ 1 i where each momenta is cal- 
culated in the laboratory frame and h + refers to either a 
K + or 7T+, is used to reject background candidates. The 
A requirement is optimized in both CS modes by max- 
imizing Ms I a N s > where Ms is the signal yield and <7/v s 
is the statistical uncertainty in Ms from the fit to the 
off-resonance data sample. The asymmetry is required 
to be less than 0.6 for both decay modes. After this fi- 
nal requirement, we find 10% and 35% improvements in 
Ms/vNs f° r CS decay modes of the D + and D+, respec- 
tively. 

Since there are differences in the mass distributions be- 
tween the data and Monte Carlo (MC) simulated @ sam- 
ples, we tune the large MC samples of generic continuum 
and BB decays, intended mainly for the accurate param- 
eterization of the peaking background under the signal. 
This background is a consequence of particle misidenti- 
fication and will be discussed in more detail later. The 
tuning procedure is as follows: the tt + {K + ) momen- 
tum scale and resolution are tuned with the D° — > ir + ir~ 
(D° — > K + K~) data sample. For the Kg momentum 
scale and resolution tuning, the D + — » if!y7r + data sam- 
ple is used. The tuning method is validated by comparing 
simulated and real data in the KgK + final state. The 
four signal decay modes are simulated and results of the 
tuning are also applied to them. 

In the branching ratio measurements, there is a peak- 
ing background due to particle misidentification. In the 
D+ — » KgK + mass region, there is a peaking struc- 
ture from D + — > Kgir + decays when a ir + is misiden- 



tified as a K + . A similar peaking structure in the 
D + — ► K { gix + mass region appears due to misidentifi- 
cation in Z?+ — * K S K + decays. The shapes and the 
yields of these peaking backgrounds are obtained from 
the tuned simulation samples and are used as the prob- 
ability density functions (PDF) for the peaking back- 
grounds. The simulated shape and normalization of the 
peaking backgrounds are checked by comparing the in- 
variant mass distributions of selected KgK + (Kgir + ) 
pairs with the K + (tt + ) mass assignment changed to a 
7r + (K + ) mass assignment. The comparison shows that 
the simulated peaking background of the tuned sample 
correctly describes these components and that misidenti- 
fication is indeed the only contribution above the struc- 
tureless combinatorial background. Uncertainties in the 
misidentification probabilities are considered as a source 
of systematic uncertainty. 

The KgK + and Kgir + invariant mass distributions af- 
ter the final selections are shown in Figs.Q]andE]together 
with the signal and background parameterizations. Clear 
signals for CF and CS decays are observed in both dis- 
tributions. The KgK + and KgiT + invariant mass dis- 
tributions are fitted using a binned maximum likelihood 
method. In all cases the signal PDF is parameterized 
using two Gaussians with a common mean value. For 
£>+ — * Kgir + , we fix the ratio of widths and the fractional 
yields in the two Gaussians because of the low statistics. 
The values of the ratio and the fraction of the broader 
Gaussian are obtained from the fit to the D + — > K s ir + 
mode and are consistent with the results of fits to MC 
simulated signal. The reduced x 2 values of the fits are 
1.8 and 2.3 for the K S K + and K s n + final states, respec- 
tively. The normalization of the mass distributions of the 
misidentified Kjix backgrounds are fixed to the values 
obtained from tuned simulation samples. Combinatorial 
background PDFs are parameterized using second and 
first-order polynomials for the K S K + and K s ir + final 
states, respectively. All the fit parameters are allowed to 
float except for the Df — > K s w + signal PDF parameters 
and the yield and the normalization of the misidentified 
backgrounds. Table U summarizes the extracted signal 
yields from the fits to data and corresponding signal ef- 
ficiencies from the simulated signal samples where final- 
state radiation has been included [10(. 



TABLE I: The extracted signal yields from the fits to data 
and corresponding signal efficiencies (e) from the simulated 
events of signal modes. The uncertainties are statistical only. 



Decay modes 


Yields 


e(%) 


D + _> K U S K+ 


100855±561 


12.59±0.01 


Df -> K%K+ 


204093±768 


13.53±0.01 


D+ -> K° s tv+ 


566285±1162 


14.19±0.01 


Dt - K° s k+ 


17583±481 


15.35±0.01 



Various contributions to the systematic uncertainties 
for the branching ratio measurements are summarized in 
Table HIl Several sources of systematic uncertainty are re- 
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M(K°K + ) (GeV/c 2 ) 



FIG. 1: Invariant mass distribution of selected K S K pairs. 
Points with error bars show the data and histograms show the 
results of the fits described in the text. Signal, D + — ► Kgir + 
background, and random combinatorial background compo- 
nents are also shown. 




M(KV) (GeV/c 2 ) 



FIG. 2: Invariant mass distribution of selected Kgic pairs. 
Points with error bars show the data and histograms show 
the results of the fits described in the text. Signal, — > 
KgK + background, and random combinatorial background 
components are also shown. The inset is an enlarged view of 
the Df region. 



duced in ratio measurements due to the similar kinemat- 
ics of CF and CS decays. Such sources include the track- 
ing and asymmetry variable efficiency differences between 
simulated data and real data. However, the systematic 
uncertainty due to particle identification efficiency does 
not cancel. The particle identification efficiency differ- 
ences between real data and simulated events are esti- 
mated independently using the decay D* + — > D°ir + fol- 
lowed by D° — » K~it + and corrections from this estimate 
are applied to signal efficiencies in Table HI These cor- 
rections are 1.000±0.007 and 0.946±0.005 for K+ and 
7r + candidates, respectively. Uncertainties in the particle 
identification corrections are included in the systematics 
estimate and are found to be 0.90% of the measured ra- 
tios. In order to validate the entire analysis procedure, 
we fit large numbers of simulated samples of generic con- 
tinuum and BB decays, and find no bias in the procedure 
within the statistical uncertainties of our measurements. 
We refit the data with various histogram binnings, differ- 
ent fit intervals, and different combinatorial background 
PDFs. We also refit the data in the D + and Df sam- 
ples separately. We estimate 0.74% and 2.00% of the 
measured ratios as the systematic uncertainties due to 
variations in fit methods for the D + and Df modes, re- 
spectively. Particle identification and the associated nor- 
malizations of the K/tt misidentified background yields 
in fits are also estimated using the measured misidenti- 
fication rates and found to be 0.16% and 0.62% of the 
measured ratio for the D + and Dj modes, respectively. 
Finally, systematic effects due to the extra constraints in 
the Df — > Kgir + signal PDF are estimated by refitting 
the data allowing the fixed parameters to change within 
their one standard deviation uncertainties. This gives a 
negligible effect in D + decays but there is a systematic 
effect corresponding to 0.37% of the measured ratio in 
.D+ decay modes. Table [TT| summarizes the systematic 
uncertainties in the branching ratio measurements. 



TABLE II: Relative systematic uncertainties in percent, 
where <r R ^ D +^ and a R i D +~. are systematic uncertainties for 

branching ratios of D + and Df decays. Sources include par- 
ticle identification (PID), fit methods, peaking background, 
and the D+ signal PDF. 



Source 


°R(D + ) (%) 


a R(D+) ( % ) 


PID 


0.90 


0.90 


Fit methods 


0.74 


2.00 


Peaking background 


0.16 


0.62 


D+ signal PDF 




0.37 


Total 


1.18 


2.31 



With the signal efficiencies and the corrections due to 
particle identification efficiency differences, we find the 
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branching ratios to be 



R(D + ) 



R(Dl 



K%K+) 



B(D+ -> K°tt+) 
0.1899 ±0.0011 ±0.0022, 
B(Df -> K° s n+) 
B(Df -> X°/f+) 
0.0803 ±0.0024 ±0.0019 



where the first uncertainties are statistical and the second 
are systematic. These are the most precise measurements 
to date and are compared to the present world average 
values in Table IIIII Our measurement of R(D + ) is in 
good agreement with previous measurements [3J and is 
larger than the naive expectation of tan 2 9c, consistent 
with the expected destructive interference effect men- 
tioned earlier. For Df decays, there is no such inter- 
ference and R(Df) is found to be greater than tan 2 9c 
by more than eight standard deviations, consistent with 
previous measurements [3j. This large deviation may be 
due to the color-suppression of the main 
amplitude. 



TABLE III: Branching ratios for the D + and the Df , and 
comparisons with previous measurements. The uncertainties 
shown combine the statistical and systematic uncertainties of 
our results. 



Branching 


Belle 


World- 


ratio 


exp. 


average [3] 


R(D + ) 


(19.0±0.2)% 


(20.6±1.4)% 


R(Df) 


( 8.0±0.3)% 


( 8.4±0.9)% 



To conclude, using 605 fb _1 of data collected with the 
Belle detector at the KEKB asymmetric-energy e + e _ col- 
lider we have measured the D + — > KgK + and Df — > 
Kgir + branching ratios with respect to the correspond- 
ing Cabibbo-favored modes. The results are B(D + — > 
K%K + )/B{D+ -> K° s it+) = 0.1899±0.0011±0.0022 and 



B(Dj 
0.0019, 



K^+)/B{Di 



K° S K+) 



0.0803±0.0024± 



the second are systematic. Using the world average 
values of CF decay rates Q, we obtain the branching 
fractions B(D+ -> K° S K+) = (2.75±0.08)x 10~ 3 and 
B(D+ -> K° s tt+) = (1.20±0.09)xl0~ 3 where the un- 
certainties are the sum in quadrature of statistical and 
systematic errors. These are consistent with the present 
world averages [H and are the most precise measurements 
to date. The ratio B(D+ -> K%K + )/B{D+ -> K° s tt+) = 
2.29±0.18 may be due to SU(3) flavor breaking and/or 
different final-state interactions in D + and Df decays. 



where the first uncertainties are statistical and 
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